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ABSTRACT 

Aims. We present an analysis of 37 high-quality extended rotation curves that highlights the existence of a new discrepancy (or a new aspect 
of an old discrepancy) between the density profiles predicted by the A Cold Dark Matter (ACDM) theory and the actual distribution of dark 
matter in galaxies. 

Methods. We compare the predicted face-value density vs. enclosed mass relationship, at large distances, to the observational data at the last 
measured radii of the rotation curves and in two whole rotation curves of high quality. A further analysis is performed by studying a relation, 
inbuilt in ACDM, that links at radius R, the enclosed halo mass MnfwW and its density p(R) in a way that is independent of the mass of the 
virialised object. 

Results. We find that the predicted density vs. enclosed mass relationship has a systematic offset with respect to the observational data. In 
test case extended rotation curves, at their last measured point, the predicted NFW densities are up to a factor 3 lower than those derived 
from the kinematics. Moreover, the abovementioned relation, inbuilt in ACDM, does not hold for the objects of our sample. Such a new outer 
discrepancy is different and maybe complementary with respect to the core/cusp issue, for which the NFW densities turn out to be higher than 
those observed and it seems to imply a global mass rearrangement of a pristine NFW-ACDM halo. 



Key words. Galaxies: kinematics and dynamics - (Cosmology:) dark matter - Galaxies: structure 



1. Introduction 

The study of disk galaxies' rotation curves has been one of 
the most successful tools for investigating the dark matter phe- 
nomenon in galaxies. Recently, the attention has been focused 
on the distribution of dark matter in galaxies, as inferred from 
the rotation curves, in particular the comparison between the 
predictions of standard A Cold Dark Matter (ACDM) the- 
ory of structure formation and observations (e.g., Salucci & 
Burkert 2000, Borriello & Salucci 2001, de Blok & Bosma 
2002, Weldrake et al. 2003, Swaters et al. 2003, Simon et al. 
2005, Gentile et al. 2004, 2005, 2006). 

More specifically, the ACDM theory predicts that dark mat- 
ter halos have a specific density distribution that follows the 
well-known NFW (Navarro, Frenk and White, 1996) profile: 



Pnfw(^) 



Ps 



(1) 



where r s and p s are the characteristic radius and density of 
the distribution. The latter is given by: 



Ps = t 



3 ln(l +c)- 



(2) 



where p c is the critical density of the Universe and A is the 
virial overdensity (see Bryan & Norman 1998). 

r s and p s are related to each other (e.g. Wechsler et al. 
2002), so eq. 1 is rather a one-parameter family of profiles, 
where the following relations link the virial mass M v ; r to the 
concentration parameter c (=r v i r /r s , where r v ; r is the virial ra- 
dius), r s and p s , at redshift z = and adapting the relations 
(similarly to Gnedin et al. 2006 and Dutton et al. 2006) to the 
cosmological parameters from the WMAP third year results 
(Spergel et al. 2006): 
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The virial radius r v i r can then be derive from r vll - — c r s . In 
ACDM, at least in a statistical sense, once the radius and the 
mass at that radius are fixed, the parameter describing the mass 
distribution (usually the virial mass M vlI ) is also known. On the 
other hand, a quite remarkable number of observations show 
that NFW profiles, displaying an inner "cusp", are inconsistent 
with data. In fact, the latter indicate profiles with a different 
characteristic, a central density "core", i.e. a region where the 
dark matter density remains approximately constant. A number 
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of studies cast doubts on the reliability of the mass modelling 
procedure and the data analysis (van den Bosch et al. 2000, 
Swaters et al. 2003, Hayashi & Navarro 2006, Valenzuela et 
al. 2007) of some galaxies. The debate is still on, despite the 
fact that most concerns are now overcome (Gentile et al. 2004, 
2005, 2006, de Blok, Bosma & McGaugh 2003, de Blok 2005). 

The simplest example of a cored halo profile is the pseudo- 
isothermal (PI, van Albada et al. 1985): 



where Rq is the core radius and po is the central density, 
that results about one order of magnitude lower than ACDM 
predictions (Donato et al. 2004, Gentile et al. 2004, 2005 and 
references therein). 

However, the above works did focus on the comparison be- 
tween ACDM and observations in the inner parts of the galax- 
ies, while it is of extreme importance to investigate also the 
outer density distribution, where ACDM halos have an equally 
strong feature: the transition between a p oc R~ l to a p oc R 
regime. In this paper, outer and inner regions have a baryonic 
perspective. The former refers to a region inside 1-2 stellar 
disk scale lengths Ru, the latter to that extended out to > 3-5Ru 
and characterised by the HI disk. Both regions, however, must 
be considered inner regions with respect to the DM distribu- 
tions. In other words, is the discrepancy between theoretical 
predictions and data present only in the inner parts of a galaxy, 
and therefore the profiles converge to the NFW one at larger 
radii, or is there a more global discrepancy extending beyond 
the NFW halo R 1 regime, as suggested also by McGaugh et 
al. (2006)? 

Note that the rotation curves at large radii, i.e. for R » 
Ru, provide a very good measure of densities and of en- 
closed dark masses. In fact, at these distances, the contribu- 
tion to the gravitational potential due to gas and stars is very 
small and the uncertainties on the accurateness/uniqueness of 
the dark-to-luminous mass decomposition play a very minor 
role: it is the rotation curve itself that essentially measures 
the physical quantities M h (R) ^ G^V 2 (R)R and p h (R) - 
l/(47rR 2 )dM h (R)/dR. Conversely, in the inner regions (R < 
Ru), usually, a complex dark-luminous mass modelling is 
needed to decompose the circular velocity into the (possibly 
equally important) dark and luminous contributions (see e.g. 
Barnes, Sellwood & Kosowsky 2004). 

2. Samples and methods of investigation 

In the present paper we investigate the outer dark matter distri- 
bution mostly resorting to two samples, heterogeneous in mass, 
covering about 3 orders of magnitude: 1) the sample of high- 
quality rotation curves selected by Donato et al. (2004), dis- 
carding the 4 galaxies with the smallest extension relative to the 
disk exponential scale length Ru and adding the galaxies: DDO 
47 (Salucci, Walter and Bordello 2003, Gentile et al. 2005) and 
ESO 287-G13 (Gentile et al. 2004); 2) a sample with rotation 
curves selected from the literature fulfilling the requirements of 
a) reaching out to at least 30 kpc, or out to 6 disk scale lengths, 
or b) with the final velocity being above 250 km s _1 . In this 



Table 1. Selected galaxies with their references. The aster- 
isk indicates that the original work provided also the dark- 
luminous decomposition of the rotation curve. 



Galaxy 


Reference 


NGC 289* 


Walsh et al., 1997 


NGC 1068 


Sofue etal., 1999 


NGC 1097 


Sofue etal., 1999 


NGC 1232* 


van Zee & Bryant, 1999 


NGC 3198* 


Blais-Ouellette et al., 2001 


NGC 3726 


Verheijen & Sancisi, 2001 


NGC 4123* 


Weiner et al., 2001 


NGC 5055 


Sofue etal., 1999 


NGC 5236 


Sofue etal., 1999 


UGC 5253 


Noordermeer et al., 2004 


NGC 5985 


Blais-Ouellette et al., 2004 


NGC 6946* 


Carignan et al., 1990 


NGC 7331* 


Bottema, 1999 


UGC 9133* 


Noordermeer et al., 2004 



way, the curves were sufficiently extended to map a region of 
the halo density profile not affected by the central slope, nor 
massively affected by the presence of the disk. It results that 
we investigate regions of halos corresponding to galactocentric 
distances out to 5% up to about 35% (the case of NGC 9133) 
of the virial radius: these zones are beyond the influence of 
the cusp, and still well into the central part of the dark matter 
halo. We selected only rotation curves that were regular out to 
the last measured radius Rf (whose average value is about 24 
kpc), obtaining a sample of 37 galaxies. In Table 1 the selected 
galaxies of sample 2 are listed, together with the corresponding 
references. 

For a number of galaxies the mass decomposition between 
the luminous and dark component of the velocity was obtained 
from the literature; we marked them with an asterisk in Table 
1. For the others, we used the method of Persic and Salucci 
(1990), see Appendix A. 

The rotation curves of sample 1 have been selected and suc- 
cessfully used (Donato et al. 2004) in order to investigate the 
core radius issue, i.e. an issue that needs more accurate data 
than the issue we want to tackle here. Since the resulting halo 
rotation curves of the galaxies of the sample are very well fitted 
by the PI profile, we will take it as representative of the dark 
matter haloes around galaxies. 

From our sample 1 we select DDO 47 and ESO 287-G1 3, as 
the best examples to show individually the outer NFW ACDM- 
data discrepancy. Then, we use the combined samples 1-2 to 
compare, at the outermost radii, the NFW ACDM predictions, 
obtained by a newly discovered structural ACDM relation, with 
the values of the mass, radius and density of the DM halos 
around galaxies. 

3. The new dark matter density discrepancy from 
the extended rotation curves of ESO 287-G13 
and DDO 47 

We analyse in detail the density profiles of the DM halos for 
the two best cases in our sample: DDO 47 (Salucci, Walter and 
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DDO 47 




Fig. 1. DM density distributions; filled circles represent the DM density inferred by the rotation curves. Solid lines corresponds 
to the NFW solution with M Nm (R f ) = M n (R f ), that give, for DDO 47, (r s , M/L B ; fl D )=(5.5 kpc, 0.5; 0.5 kpc), and for ESO 
287-G13, (r s , M/U; # D )=(17.4 kpc, 1.8; 3.3 kpc)., The dashed line corresponds to the "best fit" NFW solution, (r s , M/Li)=(19.4 
kpc, 0.7). See text for the (conservative) estimate of the uncertainties. 



Borriello 2003, Gentile et al. 2005) and ESO 287-G13 (Gentile 
et al. 2004). Note that in the present study, differently from the 
original papers, the aim is to show the existence of an outer 
observations/theory discrepancy, rather than to investigate the 
inner cusp/core issue: therefore we will analyse the data in a 
different but proper way. For consistency with the detailed anal- 
ysis made in the original papers, in this Section we used the 
relation between c and M v ; r given by Bullock et al. (2001) and 
Wechsler et al. (2002) instead of eq.[3] 

3.1. Actual dark matter halos 

We follow Salucci et al. (2003) and Gentile et al. (2004), the 
density distribution ph(R) of the dark matter halo in these galax- 
ies is given by: ph(R) = Ppi(R) (see eq. 4). Notice that in both 



objects we have also considered a Burkert halo (Burkert 1995), 
and found results coincident with the PI halo. 

3.2. NFW halos 

Let us assume that: 

M NFW (fif) = M P1 (R f ) (5) 

where Mp\{Rf) is the mass inside Rf we derive from the 
above mass model. Eq. [5] allows to derive the values for c and 
M vir : 18.4 and 6 x 10 10 M for DDO 47, and 13.3 and 7 x 10 11 
M Q for ESO 287-G13. Let us anticipate that this assumption is 
very conservative; in fact, if the mass equality expressed by eq. 
[3]occurs at radii smaller than Rf the result we claim here will be 
even more prominent; on the other hand, if the mass equality 
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DDO 47 




r/r s 

Fig. 2. Density profiles of DDO 47 assuming PNFwO^f) — 
p(Rf). See Fig.Q]for the explanation of the symbols. 

occurs at radii larger than Rf, a discrepancy from R = out 
to R — Rf will be set by definition. This assumption has also 
the desirable by-product, in the NFW framework, of implying 
reasonable values for the stellar mass-to-light (M/L) ratios. 

Given the importance of a careful analysis, we have also 
considered a different implementation of the NFW halo + disk 
+ gas mass modelling: in ESO 287-G13 we have fitted the ro- 
tation curve with this mass model by leaving M/L and M vn as 
free parameters; this latter procedure, shown in Fig. Q] yields 
very similar results with respect to those obtained by means of 
eq.0 

3.3. Test for ESO 287-G13 and DDO 47 

In Fig. Q]we compare the actual dark matter density ppi(R) and 
Pnfw(R) at the radii where the rotation curves were measured. 
A conservative estimate of the uncertainty on the density was 
derived from the mass modelling uncertainties. In fact 20% un- 
certainties in po and Rc give an average uncertainty on the den- 
sity of about 25%. A rigorous derivation goes beyond the scope 
of this paper. Then, it is clear that the density is not determined 
with the same accuracy at each radius. In fact dark matter is not 
the dominant kinematic term in the inner parts and therefore the 
density uncertainty is larger. 

We see the well-known (though debated in some case) 
cusp-core discrepancy at R — > 0, but at larger radii we now 
realise that the difference in densities changes sign and that 
the "inner" discrepancy is now reversed: in both analyses, the 
NFW halo densities, from a certain radius onward, are lower 
than the dark matter halo density. Note that such a discrep- 
ancy was also present in some of the analysis made in previ- 
ous investigations (e.g. Blais-Ouellette et al. 2001, Borriello & 
Salucci 2001, de Blok & Bosma 2002), but it was not claimed 
explicitely, neither it was investigated whether some combina- 
tion of c and M v ; r had made possible for pnfw to converge to the 
actual density inside the region mapped by the rotation curve 



data. In the present paper we instead claim that the situation 
can be described by one of these two possibilities: 1) pnfwW 
at any radius is systematically higher than the estimated density 
(so it will disagree in the core region, but will be compatible 
with the outermost density, see Fig. [2]) and the data vs. predic- 
tion discrepancy will extend to several disk scale lengths; 2) the 
discrepancy is "bivariate": actual halos (in comparison to pre- 
dicted halos) have a density deficit in the inner regions but an 
excess in the outer ones. We consider the second possibility as 
more likely since it implies stellar M/L ratios compatible with 
the predictions of stellar population synthesis models (Bell et 
al. 2003, Gentile et al. 2004). For instance, in the case of DDO 
47, imposing pnfw(^) = p{R) leads to a best-fit M/L ratio of 
0; on the other hand, with the second possibility one gets a B- 
band M/L = 0.7, within the range 0.5 - 0.8 arising from the 
observed B - V colour and the predictions of stellar population 
synthesis models (Bell et al. 2003). 

Let us stress that the availability of data at large distances is 
a crucial point: in these galaxies we just reach the radius where 
the NFW density decreases in an appreciable way. Compared 
to the discrepancy in the inner regions of galaxies, extremely 
evident and as wide as an order of magnitude, the outer dis- 
crepancy is less spectacular but not less important. 

4. The new density discrepancy from a new NFW 
halo phenomenology 

A main property of the distribution of NFW-CDM haloes 
is that M NFW (tf) = f(R, M vir , c(M vir )), while p NFW (R) - 
h(R,M vir ,c(M v i r )). From this it follows that: M NF w(R) = 
@(R,Pnfw(R), M W1I ), that is, at a fixed radius, the density and 
mass are related, though in principle in a different way in galax- 
ies of different M v ; r . The above provides us with a convenient 
way to present the structural properties of NFW halos, accord- 
ing to which observational data can be compared with the the- 
ory without knowing the virial masses of the objects. This is 
crucial because of the poor fitting performance of the NFW fits 
to rotation curves often prevents even a rough estimate of this 
quantity. The relation we find (hereafter called the M - R - p 
relation) is valid for any M wu \ this can be easily seen in Fig. [3] 
by plotting it for 3 different virial (total) masses (5 x 10 10 M , 
1 x 10 12 M Q and 1 x 10 13 M Q ) and concentration parameters 
(see eq. |2), for radii ranging from 0.001r v i r to r v i, the relation 
takes the form (see Fig . [3j : 

3 11 
log M NFW = -log pnfw(A) + —log R + 26.17 (6) 

where M NFW is in M , R in kpc and Pnfw in g cirT 3 . The 
M - R - p relation holds for more than 6 orders of magnitude 
in mass, and it is valid, for real cases, for any galactic halo at 
any radius. The scatter in c, 3 orders of magnitude in radii and 
a factor 200 in M v h introduce a negligible (< 0.05 dex) scatter 
through the 6 orders of magnitude in MnfwCR). 

The relation follows from the structural properties of NFW 
halos. It can be derived mathematically, but its physical mean- 
ing is reported to the origin of eqs. Q]and[3] Let us write the 
ratio between the enclosed mass Mnfw(< R) of an NFW halo 
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Fig. 3. The M — R — p relation: the (overlapping) solid lines refer to the ACDM haloes, for 3 different virial masses (5 x 10 10 
M , 1 x 10 12 M Q and 1 x 10 13 M ) and the 3 corresponding concentration parameters c according to Wechsler et al. (2002) (18.7, 
12.7 and 9.4, respectively). The dashed lines correspond to the + l-cr uncertainty in c taken from Wechsler et al. (2002). Radii 
range from 0.001r v ; r to r vlr . M(<R) is in M Q , p(R) in g cirT 3 and R in kpc. 



at radius r (where the density is pnfwW) and a uniform sphere 
of radius r and density pnfw(R)'- 



M NFW (< r) 
4/37rr 3 pNFw(fi) 



3(1 +cxf 



I CX \ 

ln(l+cx)-- 

V 1 + cx I 



(7) 



M NF w(<r) 



On the other hand, from eq. [6] we have: ^ — -r - 

(p NFW (R)R)- 1 ' 4 . Then, recalhng that (p NFW (R)Ry 1 / 4 = £ts££, 
one finds that the M - R - p relation is a way of expressing the 
approximation (valid for < x < 1): 



(1 +cx) l/2 ocA 



3(1 +cxf 
(c 2 x 2 ) 



(ln(l 



+ cx) 



1 + cx) 



(8) 



We compared the predicted M - R - p relation with the 
measurements of 37 galaxies. The result (Fig. |4j is that the ob- 
servational points at Rf are systematically offset compared to 
the M - R - p relation. The mean offset of ~0. 1 dex is solid 
(but note that larger offsets are also observed); in fact, the er- 
ror propagation analysis shows that a 3% uncertainty on the 
measure of the circular velocity and a 0.05 uncertainty on the 
logarithmic gradient dl °/f ^^ yield uncertainties of the order of 
0.025 dex in M and 0.06 dex in p. The uncertainties on the ob- 
jects distances are not included since they only induce random 
errors. At Rf, for a given enclosed mass M(< Rf), the observed 
density ppi{Rf) is higher than that of NFW haloes that match 
M(< Rf). The density discrepancy runs up to a factor ~3 and 
it is significant in that the adopted procedure of matching the 
CDM mass to the observed one is a conservative one. The pre- 
dicted halo profiles are based on a larger number of halos than 
the present sample of 37 galaxies, so selection or systematic 
effects cannot be completely excluded, even if the 37 galaxies 



span a large range of maximum velocities, Hubble types and 
environments. 

We tested whether the environment has an effect on our 
results: we distinguished between galaxies in "isolated" halos 
and galaxies in "subhalos". A straightforward definition does 
not exist, so we defined as "subhalos" those galaxies which ei- 
ther belong to a known group of galaxies or that have a larger 
companion within 30 R25 and 400 km s -1 . The other galaxies 
were labelled as "isolated". A more detailed investigation of 
the effect of environment goes beyond the scope of our paper, 
especially because of the heterogeneity of the data we consider 
here. Qualitatively there are no obvious differences between the 
different subsamples, even though the predicted difference be- 
tween the concentrations of halos and subhalos (Bullock et al. 
2001) cannot be excluded. 

Edge-on galaxies may represent a potential problem, as the 
Ha rotation curves could be affected by extinction (Bosma et al. 
1992) and the HI rotation curves might suffer from projection 
effects unaccounted for and from the possible lack of gas on 
the line of nodes. In Fig. [4] we have plotted the three edge-on 
galaxies (i.e., with an inclination larger than 85°) of our sample 
with a triangle; we realise that these galaxies do not lie in any 
peculiar region of the plane and they are not more discrepant 
than the other galaxies. The same holds for the galaxies from 
sample 2 for which the mass model was obtained as described 
in the previous Section. 

Baryons are expected to affect the density distribution of 
dark matter, but the effects are far from being clear: the best 
studied process is adiabatic contraction (Blumenthal et al. 
1986, Gnedin et al 2004, Sellwood & McGaugh 2005), which 
would make the halos more centrally concentrated, but mech- 
anisms with opposite effects have also been studied, such as 
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adiabatic expansion (Dutton et al. 2006) or dynamical friction 
(Tonini, Lapi & Salucci 2006). Hence, in the ACDM halos con- 
sidered here we did not take the effect of baryons into account, 
since their effect on the distribution of dark matter is still poorly 
understood. Also, baryons, when gas cools, are also expected 
to make dark matter halos more spherical (Kazantzidis et al. 
2004). In Fig.|6]we plotted the baryonic mass of the galaxies of 
our sample with mass decompositions vs. their distance from 
the relation in eq. [6] no clear correlation is observed, mean- 
ing that the discrepancy we report in the present paper is not 
straightforwardly related to the amount of baryons present in 
the galaxy, or that effects such as uncertainties in the mass dis- 
tribution create a large scatter in Fig. [6] We also note that in the 
galaxies studied here Rf » R^, i.e. Rf probes a region outside 
that most affected by adiabatic contraction. 

Warps and non-circular motions are also a potential con- 
cern for the present analysis, since in a ACDM Universe dark 
matter halos are expected to be triaxial, which would induce 
non-circular motions in the gas (Hayashi & Navarro 2006), and 
gas moving through filaments (Dekel & Birnboim 2006) will 
interact with galaxies, triggering structures on disks like those 
warps Nearly all the rotation curves collected in the present 
paper were derived using the tilted-ring fitting of the velocity 
field, which can account for warps but not for non-circular mo- 
tions. The exception is DDO 47, which was studied in detail 
by Gentile et al. (2005), using the harmonic decomposition of 
the velocity field (Wong, Blitz & Bosma 2004). So we exclude 
that warps might have an effect on the paper conclusions, while 
we cannot exlude the possibility of non-circular motions; how- 
ever, one can expect them to increase the observational scatter 
but not to have a systematic effect. 

The outer density discrepancy arises in a clear way by set- 
ting the NFW halo mass within the last point equal to the ob- 
served one at Rf\ what happens if this is valid at a different 
radius? If the radius is smaller, the density inversion happens 
at even smaller radii, and the discrepancy at Rf is worse. If 
the radius is larger, there might be no density inversion, but 
the cusp/core discrepancy (see however Section [TJ between 
ACDM predictions and observations would be increased and 
would be present at any radius where baryons are present. The 
latter case is shown in Fig. [2] where at the last point we im- 
posed ph(Rf) = Pnfw(R{) instead of Mh(Rf) = M^^iRf): the 
inner discrepancy becomes worse. Even though the method of 
analysis is different, this effect might be related to the results 
of Seigar et al. (2006), whose NFW models (including possi- 
ble adiabatic contraction) match the observed rotation curves 
of two galaxies in their outer parts (~ 10 kpc), but overesti- 
mated the inner rotation curves. One of the two cases shown 
by Seigar et al. is a barred galaxy and its fitted concentration 
parameter is small compared to the averaged predicted value; 
the authors discuss possible evidence for the absence of adia- 
batic contraction or alternatively of some dynamical effects that 
compensates for adiabatic contraction. 

By assuming the Burkert halo instead of the NFW halo 
leads to the results shown in Fig. [5] Here we consider the 
one-parameter family of halos defined by the Burkert halo 
and the empirical relation between the central density po and 
the core radius rrj found by Salucci & Burkert (2000): po = 



3 x 10~ 24 (ro/kpc)~ 2 / 3 g cirT 3 . The virial masses were computed 
by integrating the density profile until the mean density was 
A times p c , where A was derived following Bryan & Norman 
(1998). In this case (see Fig.[5]l we find that halos with differ- 
ent masses do not overlap (because for the chosen axes there 
is no such an approximation as eq. [8]), the agreement with the 
observations is much better than in the case of the NFW halos. 
This means that the Burkert halo, known to fit the inner parts 
rotation curves better, is also a better representation of the ob- 
servations at the last radii probed by rotation curves (~ 5-35% 
of the virial radius). 

Let us point out that while the present work shows that a 
theory vs observations discrepancy extends well beyond the 
very inner regions of spirals, it leaves open the possibility that, 
for radii r > 0.2 - 0.3r v i r , the density of the DM around galax- 
ies converges to a NFW profile. This was hinted at by Prada 
et al. (2003) who have investigated the very outer dark matter 
profile by means of the kinematics of satellites around isolated 
galaxies. However, they probe an outer region with respect to 
those considered here and in a very low resolution mode: their 
projected radii range from 20 to 350 kpc, in 100 kpc bins. 
Brainerd (2004a,b), using weak galaxy lensing in addition to 
the dynamics of satellite galaxies, reach similar conclusions 
for scales > 50 h~ l kpc. The results showed in the present pa- 
per hold for galaxies but the situation on the scales of galaxy 
clusters might be different: indeed, Vikhlinin et al. (2006) and 
Zappacosta et al. (2006) show agreement between X-ray data 
and ACDM mass profiles. 

5. Conclusions 

An accurate mass modelling of the external regions in the case 
of a couple of test-case spirals and a careful determination of 
the densities and enclosed masses of the dark matter haloes 
at the farthest radii at which 37 high quality rotation curves 
have been measured, has brought to the discovery of a new 
problem/discrepancy for the ACDM/NFW haloes. 

In fact, in addition to the well-known evidence for which 
in the inner regions of galaxies (R < 2Ro) the DM haloes 
show a flatfish density profile, with amplitudes up to one or- 
der of magnitude lower than the ACDM predictions, at outer 
radii (R > 4Rd) the measured DM halo densities are found 
higher than the corresponding ACDM ones. This implies ei- 
ther that the shallow-steep disagreement extends all over 1-2 
times the galaxy optical radius (which would be the case if 
PNFw(^f) = p(R{)) or that there is a complex data vs. theory 
disagreement. 

While the statistical significance and the level of the dis- 
crepancy must be investigated with more and outer data, there 
is already an evidence for this discrepancy in most galaxies 
with high quality data. The DM halo density, known to have a 
core in the internal regions, does not seem to converge to the 
NFW profile at 4-6 R D . This implies an issue for ACDM that 
should be investigated in future, when, due to improved obser- 
vational techniques, the kinematic information will be extended 
to the ~ 100 kpc scale (Gentile et al., in prep.). 

This new discrepancy provides additional information on 
the nature of the cusp/core issue: self-interacting or annihilat- 
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Log(p 3/4 (R) R 11/4 ) 

Fig. 4. The same as Fig [3] zoomed to the radial range from 1% to 20% of the virial radius. Filled circles are the observational 
data of samples 1 and 2 at Rf. The three edge-on galaxies are denoted by empty circles. Empty triangles denote the objects of 
sample 2 without published mass modelling. Red symbols are isolated halos and black symbols are subhalos, according to our 
definition in the text. 




-17 -16 -15 -14 

Log(p 3/4 (R) R 11/4 ) 

Fig. 5. The same as FiglH but showing that Burkert halos (Salucci & Burkert 2000) predict a mass dependent Mvs(p a R b ) 
relation in good agreement with observations. Symbols are the same as Fig |H and the 3 lines correspond to 3 different virial 
masses (5 x 10 1() M , 1 x 10 12 M and 1 x 10 13 M Q , see text). 



ing dark matter proposed as the cause for the inner discrepancy 
may be in difficulties in that it will cause a rapid convergence 
to the NFW profile in the luminous parts of galaxies and be- 
yond once a critical density value is reached. The discrepancy 
points to a scenario of modified CDM profiles, to a global mass 
or angular momentum rearrangement (e.g. Dutton et al. 2006, 
Tonini et al. 2006) that would remove dark matter from the in- 
nermost parts to the radii probed by the outermost regions of 
rotation curves. 



References 

Barnes, E.I., Sellwood, J.A., Kosowsky, A., 2004, AJ, 128, 2724 
Bell, E.F., Mcintosh, D.H., Katz, N., Weinberg, M.D., 2003, ApJS, 
149, 289 

Blais-Ouellette, S., Amram, P., Carignan, C, 2001, AJ, 1952, 1964 
Blais-Ouellette, S., Amram, P., Carignan, C, Swaters, R., 2004, A&A, 
420, 147 

Brainerd, T. G., 2004a, preprint (astro-ph/0409381) 
Brainerd, T. G., 2004b, AIP Conf. Proc. 743: The New Cosmology: 
Conference on Strings and Cosmology, 743, 129 



8 



Gianfranco Gentile, Chiara Tonini and Paolo Salucci: Outer ACDM Halo Density Profiles 



"1 0.4 

CD 
+ 

I 0.2 



O 



t — i — i — I — i — i — r 



• : 



_i i i i i i i i_ 



10 



i i 



12 



M 



baryons 



Fig. 6. Distance from eq. [6] vs. baryonic mass of the galax- 
ies with mass decompositions: there is no clear correlation be- 
tween the two quantities. Units are the same as Fig. [3] 



Blumenthal, G. R., Faber, S. M., Flores, R., & Primack, J. R., 1986, 

ApJ, 301, 27 
Bordello, A., Salucci, P., 2001, MNRAS, 323, 285 
Bosma, A., Byun, Y, Freeman, K. C, Athanassoula, E., 1992, ApJ, 

400, L21 
Bottema, R., 1999, A&A, 348, 77 
Bryan, G.-L., Norman, M.-L., 1998, ApJ, 495, 80 
Burkert, A., 1995, ApJ, 447, L25 

Carignan, C, Charbonneau, P., Boulanger, F., Viallefond, F, 1990, 

A&A, 234, 43 
de Blok, W. J. G., Bosma, A., 2002, A&A, 385, 816 
de Blok, W. J. G., Bosma, A., McGaugh, S. 2003, MNRAS, 340, 657 
de Blok, W. J. G. 2005, ApJ, 634, 227 
Dekel, A., Birnboim, Y., 2006, MNRAS, 368, 2 
Donato, F, Gentile, G, Salucci, P., 2004, MNRAS, 353, L17 
Dutton, A. A., Courteau, S., de Jong, R., Carignan, C, 2005, 619, 218 
Dutton, A. A., van den Bosch, F. C, Dekel, A., Courteau, S., 2006, 

preprint ( |astro-ph/0604553| l 
Gentile, G., Salucci, P., Klein, U., Vergani, D., Kalberla, P., 2004, 

MNRAS, 351, 903 
Gentile, G, Burkert, A., Salucci, P., Klein, U., Walter, F, 2005, ApJ, 

634, L145 

Gentile, G, Salucci, P., Klein, U., Granato, G. L., 2006, MNRAS in 

press, astro-ph/061 1355 
Giraud, E., 1998, AJ, 116,2177 

Gnedin, O. Y, Kravtsov, A. V., Klypin, A. A., & Nagai, D., 2004, ApJ, 
616, 16 

Gnedin, O. Y, Weinberg, D. H., Pizagno, J., Prada, F, Rix, H.-W., 

2006, preprint ( astro-ph/0607394 1 
Hayashi, E., Navarro, J. F, 2006, MNRAS, in press, |astro-ph/0608376| 
Kazantzidis, S., Kravtsov, A. V., Zentner, A. R., Allgood, B., Nagai, 

D., Moore, B., 2004, ApJ, 611, L73 
McGaugh, S. S., de Blok, W. J. G, Schombert, J. M., Kuzio de Naray, 

R., & Kim, J. H. 2006, ApJ, in press Jastro-ph/06 12410) 
Navarro, J.F, Frenk, C.S., White, S.D.M., 1996, ApJ, 462, 563 
Noordermeer, E., van der Hulst, T., Swaters, R., 2004, Proceedings 

of "Baryons in Dark Matter Halos". Novigrad, Croatia, 5-9 Oct 

2004. Eds: R. Dettmar, U. Klein, P. Salucci. Published by SISSA, 

Proceedings of Science, http://pos.sissa.it, p. 68. 
Persic, M., Salucci, P., 1990, MNRAS, 247, 349 



Prada, F, Vitvitska, M., Klypin, A., Holtzman, J. A., Schlegel, D. J., 
Grebel, E. K, Rix, H.-W., Brinkmann, J., McKay, T. A., Csabai, 
I., 2003, ApJ, 598, 260 
Salucci, P., Burkert, A., 2000, ApJ, 537, L9 
Salucci, P., Walter, F, Borriello, A., 2003, A&A, 409, 53 
Seigar, M. S., Bullock, J. S., Barth, A. J., & Ho, L. C. 2006, ApJ, 645, 
1012 

Sellwood, J. A., McGaugh, S. S. 2005, ApJ, 634, 70 
Simon, J.D., Bolatto, A.D., Leroy, A., Blitz, L., Gates, E.L., 2005, 
ApJ, 621,757 

Sofue, Y, Tutui, Y, Honma, M., Tomita, A., Takamiya, T., Koda, J., 

Takeda, Y, 1999, ApJ, 523, 136 
Spergel, D. N., et al., 2006, preprint jastro-ph /0603449 ) 
Swaters, R.A., Madore, B.F, van den Bosch, EC, Balcells, M., 2003, 

ApJ, 583, 732 
Tonini, C, Lapi, A., Salucci, P., 2006, ApJ, 649, 591 
Valenzuela, O., Rhee, G, Klypin, A., Governato, F, Stinson, G., 

Quinn, T., Wadsley, J. 2007, ApJ, in press, |astro-ph/0509644 1 
van Albada, T.S., Bahcall, J.N., Begeman, K, Sancisi, R., 1985, ApJ, 

295, 305 

van den Bosch, EC, Robertson, B.E., Dalcanton, J., de Blok, W.J.G., 

2000, AJ, 119, 1579 
van Zee, L., Bryant, J., 1999, AJ, 118, 2172 
Verheijen, M. A. W., Sancisi, R., 2001, A&A, 370, 765 
Walsh, W., Staveley-Smith, L., Oosterloo, T., 1997, AJ, 113, 1591 
Wechsler, R.H., Bullock, J.S., Primack, J.R., Kravtsov, A.V., Dekel, 

A., 2002, ApJ, 568, 52 
Weiner, Benjamin J., Williams, T. B., van Gorkom, J. H, Sellwood, J. 

A., 2001, ApJ, 546,916 
Weldrake, D.T.F, de Blok, W.J.G., Walter, F, 2003, MNRAS, 340, 12 
Wong, T., Blitz, L., Bosma, A., 2004, ApJ, 605, 183 

Appendix A: Persic & Salucci (1990) method to 
derive the disk mass 

The total velocity profile is given by assuming the total cen- 
trifugal equilibrium between the two components, the baryons 
(only the disk is taken into account, as the bulge just affects 
the centre of the rotation curve (RC), and the stellar halo mass 
and the HI disk are negligible for our purposes) and the Dark 
Matter: 



V\R) = Vl{R) + VliR) 



(A.l) 



The halo and disk mass inside a given radius R is Mu(R) = 
G 1 V^(R) R and the disk contribution to the circular ve- 
locity is: V^(R) = GMuf/R respectively, where / = 
1/2 (R/R D ) 2 (IqKq - I\Ki)\ m R D is defined in terms of the 
modified Bessel functions /„ and K„. Let us now consider 
the moments of equation ( IA.U . and let us define the log- 
arithmic slopes of the total, disk and halo velocity respec- 



tively: V = 



d log V(R) 
d log R 



V H = 



d log V H (R) 
d log R 



V = 
where R D?t 



d log V D (R) 
d log R 



-0.27, 



3.2R D . While V is ob- 



served, Vh is unknown, and it is related to the Dark Matter 

2V H + 1. Since 



distribution by: d(log M H {R))/d{log R) Ropl 

tow write 

ratio at the optical radius R opt : 



Md ~ l.lVp(R)R/G. we can now write the disc-to-total mass 



V H -V 



M D (R opt ) 



M total (fi opt ) Vh+O.IIV + 0.30 



(A.2) 
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It is reasonable to assume that at the optical radius the halo 
RC is not decreasing, and the DM density is not increasing, 
therefore < V H < 1 ; in the maximum disk hypothesis, we 
can find the value of Vh in the interval [0, 1] that maximises 
the ratio of equation CO: 4Md( ^ Mm "' (J? " pi)1 = and in this 



case it is easy to see that Vh = 1 . From equation jA.2j we then 
obtain: 

M D - 1.32 V 2 opt R opt (0.77 - 0.83V) (A.3) 

which allows us to obtain the halo velocity profile at any given 
point in its luminous and dark contributions, by means of Md, 
/ and eq. IA.1I 

Notice that this minimum halo assumption used to derive 
the disk mass might not be completely correct and it might en- 
hance the evidence for an inner core, but it is perfectly legiti- 
mate in the present procedure in that in any case it just (slightly) 
under-estimates the DM density at large radii, bringing even 
more support to our claim. 

In any case the disk mass, obtained by means of eq. IA.3I 
plays a very minor role in the results of this work, in fact at 
large distances the mass and density of the dark matter at any 
outermost radius (denoted as Rf) do not depend significantly on 
this quantity. We have, supposing to deal with a system in virial 
equilibrium: 



Pft) = , —5 1+2 



4nGR\ \ dlogR 



+ ^-f(R f /R D ) (A.4) 

Rr, 



(see Fall & Efstathiou 1980). We realise that, at most, the sec- 
ond term is a correction of 30%. So errors (even quite large) in 
the estimate of M D little propagate in the estimate of p(Rf). 



